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Abstract 

Considering the diffraction of a plane wave by a periodically corrugated half-space, we show 
that the transformation of the refracting medium from positive/negative phase-velocity to neg- 
ative/positive phase-velocity type has an influence on the diffraction efficiencies. This effect in- 
creases with increasing corrugation depth, owing to the presence of evanescent waves in the troughs 
of the corrugated interface. 
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I. INTRODUCTION 

The rediscovery of isotropic dielectric-magnetic materials exhibiting phase velocity 
vector opposed in direction to the time-averaged Poynting vector has prompted a flurry of 
publications during the last four years [2 L 3 L Many interesting effects have been predicted, 
with some experimental backing as well 4 , g, (| . 

Though several names have been proposed for this class of materials, we think that the 
most descriptive is: negative phase-velocity (NPV) materials 7]. In contrast, the phase 
velocity and the time-averaged Poynting vectors are co-parallel in positive phase-velocity 
(PPV) materials. PPV materials are, of course, commonplace and require no introduction. 

That the intrinsic difference between NPV and PPV materials has recognizable conse- 
quences is easily gauged from a simple problem: reflection and refraction of a plane wave due 
to a homogeneously filled half-space. Let vacuum be the medium of incidence, while e 2 and 
fi2 denote the relative permittivity and relative permeability of the medium of refraction. 
Let a linearly plane wave be incident on the planar interface of the two mediums at an angle 
$0) (l^ol < 7r /2), from the normal to the interface, and p(6o) be the reflection coefficient. If 
the transformation {e 2 — > — e^, P2 —* — /4i is implemented, then p(9 ) — > p*(0 ), where the 
asterisk denotes the complex conjugate [8||. Thus, the replacement of a NPV/PPV half- 
space by an analogous PPV/NPV half-space changes the phase of the reflection coefficient 
but not its magnitude. 

What wou.d happen if the iaterfaee were to he eclated Q 7 Surfaced gr ati ngS 
are periodically corrugated surfaces that are commonly used in electromagnetics, and many 
theoretical techniques are available to compute their diffraction efficiencies [10]. Therefore, 
we decided to compute and compare the diffraction efficiencies of PPV and NPV surface- 
relief gratings. In this report, we present our chief results here. Section II contains a sketch 
of the theoretical method we chose, while Section III is a discussion of the numerical results 
obtained. An exp(— iut) time-dependence is implicit. 

II. THEORY 

In a rectangular coordinate system (x,y,z), we consider the periodically corrugated 
boundary y = g(x) = g(x + d) between vacuum and a homogeneous, isotropic, linear 
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material, with d being the corrugation period. The region y > g(x) is vacuous, whereas 
the medium occupying the region y < g(x) is characterized by complex-valued scalars 
£2 = £2R + i^2i an d A*2 = f^2R + ^21- If this medium is of the NPV type, then [?], 11 1 



E2i?M + A^M < 0; (1) 

otherwise 

e2i?M + A*2Jj|e 2 | > . (2) 

A linearly polarized electromagnetic plane wave is incident on this boundary from the region 
y > g(x) at an angle 9 , (\9 \ < n/2), with respect to the y axis. 

Let the function f(x,y) represent the z-directed component of the total electric field for 
the s-polarization case, and the z-directed component of the total magnetic field for the 

n 

p-polarization case [12j . Outside the corrugations, f(x,y) is rigorously represented by the 
following Rayleigh expansions 

f(x, y) = exp \i (a x - f3^y) + 

^2 Pn exp [i (a n x + fifty)] , y > maxg(x) , (3) 



ffav) = ^2 T n exp [i (a n x - fifty)] , y < mmg(x) . (4) 

n=— oo 

Here, {pn} n= t^ and {7n} n= _rc are scalar coefficients to be determined; and 

a = ^ sin 6>o 
a n = a + 2nii/d 

= ' <5> 

where c is the speed of light in vacuum. Note that fift is either purely real or purely 
imaginary; and the conditions 
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are appropriate for plane waves in the vacuous half-space y > g(x). The refracting half- 
space y < g(x) being filled by a material medium, e 2 / > and ^21 > by virtue of causality. 
The refracted plane waves must attenuate as y — > —00, which requirement leads to the 
condition 

Im [/#>] > . (7) 



Fulfilment of this condition automatically fixes the sign of Re 



/ft 



(2) 



regardless of the signs 

of e 2 R and /z 2 r- We must note here that the transformation {e 2 _R — > —e 2 R, H2R ~^ ~^2r} 

(2) 

alters the signs of the real parts of all flk . 

Boundary conditions at y = g(x) require the continuity of the tangential components of 
the total electric field and the total magnetic field. Hence, 



f(x,g(x)+) = f(x,g(x)-) 

n ■ V/(x, g(x)+) =a- l h- V/(x, g(x)-) 



(8) 



where a = fi 2 for the s-polarization case and a = e 2 for the p-polarization case, while n is 
a unit vector normal to the boundary. 

At this stage we invoke the Rayleigh hypothesis [^| — that is, we assume that expansions 
(J2J) and (jlj), which are strictly valid outside the corrugated region, can be used in the bound- 
ary conditions (JBJ). Doing so, and then projecting into the Rayleigh basis {exp(i oc m x)} m *^, 
we obtain a system of linear equations for {p n } n= ^ and {Tn} n= ^,- Following Maradudin 
[13! p. 427], we write down the system in matrix form as 



(9) 



Mn M12 




K 




U 


M21 M22 




T 




V 



The (m, n)-th elements of the four matrixes on the right side of © are 



Mn 
M12 
M11 
M22 



~D mn (f3n , 



D mn ( f3n ) 



/3 n ^ D mn (/3n ^ ) OL n E mn ^ ) 

— f3 n D mn ( f3 n ) 4" Ot n E mn ( (5 n , 



(10) 
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while the m— th elements of the four column vectors in the same equation are 



R 
T 
U 
V 



D mQ {-f3 { Q l) ) + a E m0 {-f3 { l) ) 



> . 



(11) 



The integrals appearing in the foregoing equations are defined as 

cd 



■Dm.n. ( U J 



and 



E mn (u) 



— J exp [— i — (m — n) x + iug(x)\ dx 

1 2,71 

— \ g\x) exp [— i — (m — n) x + iug(x)] dx . 
d Jo d 



(12) 



(13) 



with the prime denoting differentiation with respect to argument. 

Equation (JSJ) has to be appropriately truncated and solved to determine the reflection 
coefficients p n and refraction coefficients r„. Diffraction efficiencies 



Re 



0; 



(1) 



(i) 



\Pn\ 2 , 



(14) 



are defined for the reflected orders. The normalized power absorbed across one period of 
the corrugated interface is given by 



Re 



Po a n,m 



Pn ^ D mn ( J P^ ^ 



+ 



(15) 



The principle of conservation of energy requires that 



(16) 



When we implemented the procedure presented, we checked that the condition (fTB*j) was 
satisfied to an error of 10 ppm. This was usually achieved by retaining 15 terms (i.e., 
—7 < n < 7) in the Rayleigh expansions and (j3J) of the fields. 
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III. NUMERICAL RESULTS AND DISCUSSION 



We chose the corrugations to be sinusoidal: g(x) = 0.5 h cos(27rx / d) . For this type 
of boundary between vacuum and a penetrable dielectric medium, good results have been 
obtained for h/d < 0.3 14, We calculated diffraction efficiencies for refracting mediums 
of both the PPV (e 2 = 5 + iO.01, // 2 = 1 + iO.Ol) and the NPV (e 2 = -5 + iO.01, /x 2 = 
— 1 + iO.01) types. Calculations were made for both the s- and the p-polarization cases. 
Fixing the ratio uod/c = 27r/l.l, we plotted the diffraction efficiencies e r and e r _ x as well as 
the absorption P a as functions of 9 G [0, 7r/2) for h/d = 0.07 (Figure [IJ, h/d = 0.14 (Figure 
I2J) and /i/d = 0.21 (Figure ED- 

When h/d = — i.e., when the interface is planar — it has been shown 8] that the 
transformation {e 2 # — > —^2r,^2R —^2r} does not change eg. No wonder, the same trans- 
formation does not seem to be very effective in affecting e r Q when h/d = 0.07. As the 
corrugations grow deeper (i.e., as h/d increases in value), the presented data shows that the 
transformation of the refracting medium from NPV/PPV to PPV/NPV increasingly affects 
e r Q and P a . 

Why should this be so? Now, for a planar interface, the transformation 
{e 2 R — > — t2R, l^2R — ► — I^2r\ leaves the magnitude of the reflection coefficient only unchanged 
for non-evanescent incident plane waves; but that is not a true statement for incident evanes- 
cent plane waves 1(|. In the troughs of the corrugated interface, the total field that exists 
has both specular (n = 0) and nonspecular (n 7^ 0) components. Most of the nonspecular 
components are like evanescent plane waves because they are characterized by Re (5n^ = 0. 
Their presence ensures that the diffraction efficiencies are affected by the transformation of 
the refracting medium from NPV/PPV to PPV/NPV. 

Before concluding, let us point out that the numerical results presented here for NPV 
surface-relief gratings agree with the results of a perturbational approach, thereby validating 
the limited use of the Rayleigh hypothesis for NPV gratings in the same way as for PPV 
gratings Q|. Also, the emergence of homogeneous NPV materials promises new types of 
gratings which could be significantly different from their PPV counterparts. 
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angle of incidence (deg) 



FIG. 1: Diffraction efficiencies e r and e r _ 1 as well as the normalized absorbed power P a as functions 
of the incidence angle 9q, for a sinusoidally corrugated interface between vacuum and a linear 
homogeneous medium. The interface function g(x) = 0.5 hcos(2irx/d), where h/d = 0.07 and 
Lodjc = 27r/l.l. The refracting medium is of either the PPV (e2 = 5 + i0.01, ^2 = 1 + 20.01) or the 
NPV (e2 = — 5 + iO.Ol, fj^ = —1 + iO.Ol) type. Calculations were made for both the s- and the 
p-polarization cases. 
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FIG. 2: Same as Figured but for h/d = 0.14. 
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